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Dihadron spectra in high-energy heavy-ion collisions are studied within the next-to-leading order 
perturbative QCD parton model with modified jet fragmentation functions due to jet quenching. 
High-pT back-to-back dihadrons are found to originate mainly from jet pairs produced close and 
tangential to the surface of the dense matter. However, a substantial fraction also comes from jets 
produced at the center with finite energy loss. Consequently, high-pT dihadron spectra are found to 
be more sensitive to the initial gluon density than the single hadron spectra that are more dominated 
by surface emission. A simultaneous ^-Ht to both the single and dihadron spectra can be achieved 
within a range of the energy loss parameter eo = 1.6 — 2.1 GeV/fm. Because of the flattening of the 
initial jet production spectra at y/s = 5.5 TeV, high pr dihadrons are found to be more robust as 
probes of the dense medium. 
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Jet quenching as discovered in high-energy heavy-ion 
collisions at the Relativistic Heavy-ion Collider (RHIC) 
is manifested in both the suppression of single inclusive 
hadron spectra at large transverse momentum (px) [J 
and the disappearance of the typical back-to-back jet 
structure in dihadron correlation in vacuum Since 
jet quenching is caused by parton energy loss which in 
turn depends on the gluon density and transport coeffi- 
cient of the medium, detailed study of the suppression of 
large px hadron spectra and correlations can shed light 
on the properties of the dense medium 0] . 

In heavy-ion collisions, the spatial distribution of the 
initial jet production points is given by the nuclear over- 
lap, T AB (b,r) = t A (r)t B (\b - r|) with t A (r) being the 
thickness function of each nucleus. The suppression fac- 
tor of the leading hadrons from jet fragmentation will 
depend on the total parton energy loss which in turn is 
related to the weighted gluon density integrated along the 
jet propagation path. Therefore, measurements of large 
Pt hadron suppression can be directly related to the av- 
eraged gluon density and medium transport coefficient. 
As the average gluon density increases with colliding en- 
ergy and centrality, one should expect continued decrease 
of the suppression factor until particle production in the 
outer corona of the medium dominates the single hadron 
spectra In this case, the suppression factor of high- 
Pt hadron spectra will lose its effectiveness as a sensitive 
probe since it only depends on the thickness of the outer 
corona which varies very slowly with the initial gluon 
density. In this Letter we will employ for the first time a 
next-to-leading order (NLO) perturbative QCD (pQCD) 
parton model to study the suppression of both single and 
dihadron spectra due to jet quenching. In particular, we 
will investigate the robustness of back-to-back dihadron 
spectra as a probe of the initial gluon density when single 
hadron spectra suppression become fragile. 

Within the NLO pQCD parton model, large px hadron 
production cross section in p + p collisions can be ex- 



pressed as a convolution of NLO parton-parton scattering 
cross sections, parton distributions inside nucleons and 
parton fragmentation functions (FF). The calculations 
discussed in this Letter are carried out within a NLO 
Monte Carlo based program Q which utilizes two-cutoff 
parameters, 5 S and S c , to isolate the soft and collincar 
divergences in the squared matrix elements of the 2^3 
processes. The regions containing the divergences are 
integrated over in n-dimensions and the results are com- 
bined with the squared matrix elements for the 2 —> 2 
processes. This results in a set of two-body and three- 
body weights, each of which depends on the cut-offs. 
However, this dependence cancels when the weights are 
combined in the calculation of physical observables. 

For the study of large px single and dihadron pro- 
duction in A + A collisions, we assume that the initial 
hard scattering cross sections are factorized as in p + p 
collisions. As in the lowest-order (LO) pQCD parton 
model study we further assume that the effect of 
final-state interaction between produced parton and the 
bulk medium can be described by the effective medium- 
modified FF's, 
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where z' c = pt/(pt c - A£ c ), z' g = (L/X)p T /AE c are 
the rescaled momentum fractions, AE C is the average 
radiative parton energy loss and (L/X) is the number of 
scatterings. The FF's in vacuum Z?^ c (z c ,^ 2 ) are given 
by the KKP parameterization 0]. 

According to recent theoretical studies [1, [l(| the 
total parton energy loss in a finite and expanding medium 
can be approximated as a path integral, 
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for a parton produced at a transverse position r and trav- 
eling along the direction n. (dE/dL)^ is the average 
parton energy loss per unit length in a 1-d expanding 
medium with an initial uniform gluon density po at a 
formation time tq for the medium gluons. The energy 
dependence of the energy loss is parameterized as 

dE 
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from the numerical results in Ref. Q in which thermal 
gluon absorption is also taken into account. The param- 
eter eq should be proportional to po- 

Neglecting transverse expansion, the gluon density dis- 
tribution in a 1-d expanding medium in A + A collisions 
at impact-parameter b is assumed to be proportional to 
the transverse profile of participant nucleons , 

p g (r, b, r) = Z^^&[ tA (p) + t A (\b - r|)]. (4) 

The average number of scatterings along the parton prop- 
agating path is 
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Here, we neglect the time-dependence of the cross sec- 
tion and characterize it by the mean free path Ao via 
do = l/(poAo) and a hard-sphere nuclear distribution is 
used. The parton distributions per nucleoli f a / A (x a ,ic) 
inside the nucleus are assumed to be factorizable into 
the parton distributions in a nuclcon given by CTEQ6M 
parameterization [TTJ] and the new HIJING parameter- 
ization [l2| of the impact-parameter dependent nuclear 
modification factor, including the isospin dependence. 

The calculated single inclusive 7r° and ir + +ir~ spectra 
in the NLO pQCD parton model in p + p collisions agree 
well with the experimental data at the RHIC energy with 
the factorization scale in the range \i = 0.9 ~ 1.5pr [H[. 
We use the same factorization scale in both p + p and 
A + A collisions in our calculation. Shown in Fig. [1] are 
the nuclear modification factors, 
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for single inclusive tt° spectra calculated in both leading- 
order (LO) and NLO pQCD parton model as compared 
to the PHENIX data on central Au + Au collisions at 
= 200 GeV. The factorization scale in the NLO result 
is set (i = l.lpT though the nuclear modification factor 
R AA is not at all sensitive to the choice of [i. However, 
R AA in NLO calculation is always smaller than the LO 
result because of the relative larger ratio of gluon/quark 
jets in NLO than in LO calculation and gluon energy loss 
is assumed to be 9/4 larger than that of a quark. In both 
calculations, we have chosen the parameters as /iq = 1.5 
GeV, eo^o = 0.5 GeV and To = 0.2 fm/c. The results are 
not sensitive to small values of tq . 
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FIG. 1: (color online). Nuclear modification factors for 7r° 
spectra in LO (dashed) and NLO pQCD (solid) in Au+Au(0— 
10%) collisions at y^i = 200 GeV with different values of 
energy loss parameter eo = 1.48, 1.68 and 2.08 GeV/fm (from 
top to bottom) compared with data 



Because of jet quenching, the dominant contribution 
to the measured single hadron spectra at large pr comes 
from those jets that are initially produced in the outer 
corona of the overlap region toward the direction of the 
detected hadron. This is clearly illustrated in Fig. [5] by 
the spatial distribution of the production points of those 
jets that have survived the interaction with the medium 
and contribute to the measured spectra. Jets produced 
in the region away from the detected hadron are severely 
suppressed due to their large energy loss and don't con- 
tribute much to the final hadron spectra. As pointed out 
in Ref. Q, when jets produced in the inner part of the 
overlapped region are completely suppressed due to large 
initial gluon density, the final large px hadron produc- 
tion is dominated by "surface emission" . High-p^ hadron 
yield via such surface emission should be proportional to 
the thickness of the outer corona which decreases with 
the initial gluon density. Therefore, the suppression fac- 
tor for single hadron spectra should continue to decrease 
with the initial gluon density as shown in Fig. [3J The 
variation, nevertheless, is very weak when surface emis- 
sion becomes dominant and single hadron suppression is 
no longer a sensitive probe of the initial gluon density. 

To find robust probes of the high initial gluon density 
when single hadron spectra suppression becomes fragile, 
we will study back-to-back dihadron spectra within NLO 
pQCD parton model in this Letter. To quantify dihadron 
spectra, a hadron-triggered fragmentation function, 
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was introduced la as a function of z T =pf s °/ p " ls , which 
is essentially the away-side hadron spectrum associated 
with a triggered hadron. Both hadrons are limited to the 
central rapidity region |y tri s> asso | < 0.5 and the azimuthal 
angle relative to the triggered hadron is integrated over 
|A0| > 2.5. 
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FIG. 2: (color online). Spatial transverse distribution (arbi- 
trary normalization) of the initial parton production points 
that contribute to the final 7r° at 8 < pr < 15 GeV along 
<f> — 7r/2. The insert is the same distribution projected onto 
the y-axis. 



As in the NLO calculation of single hadron spectra, 
one also has to fix the factorization scale in the NLO 
calculation of dihadron spectra, which is chosen to be 
the invariant mass of the dihadron M 2 = (pi+p2) 2 - The 
NLO results on associated hadron spectra D pp (zt,Pt 16 ) 
in p+p collisions are compared to the d+Au data at y/s = 
200 GeV in Fig. [4] assuming no nuclear effects in d + Au 
collisions. The overall normalization of NLO results is 
quite sensitive to the factorization scale as indicated by 
the shaded region corresponding to fx = 0.8 — 1.8M. We 
will use /i = 1.2M in this study with which NLO results 
describe the experimental data well. We use the same 
scale in the calculation of dihadron spectra in Au + Au 
collisions and the results for Daa(zt,p^ 1s ) agree well 
with STAR data as shown in Fig.[4]using the same energy 
loss parameter eo = 1-68 — 2.08 GeV/fm. The centrality 
dependence of the data is also reproduced by the NLO 
calculation [l3j]. The nuclear modification factor of the 
hadron-triggered fragmentation function 
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as plotted in the lower panel is coincidcntally similar to 
the modification factor for single hadron spectra Raa- 

We have adjusted the energy loss parameter eo to fit 
experimental data on both single (do aa/ 1 dprdy between 
pr = 4 — 20 GeV/c) and away-side spectra \Daa{zt)] in 
the most central Au + Au collisions at s/s = 200 GeV. 
The best fits occur for e = 1-6 — 2.1 GeV/fm as shown by 
the x 2 in Fig- El The fact that both x 2 's reach their min- 
ima in the same region for two different measurements is 
highly nontrivial, providing convincing evidence for the 
jet quenching description. 
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FIG. 3: (color online). The suppression factors for single 
(Raa), dihadron (Iaa) spectra at fixed transverse momen- 
tum and x 2 /d.o.f. (curves with filled squares) in fitting ex- 
perimental data on single [l4| (pr = 4 — 20 GeV/c) and away- 
side spectra [3 (p^ = 8-15 GeV, z T = 0.45 - 0.95) in 
central Au + Au collisions at y^s = 200 GeV as functions of 
the initial energy loss parameter eo- 
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FIG. 4: (color online). Hadron-triggered fragmentation func- 
tions Daa(zt) and the medium modification factors Iaa(zt) 
in NLO pQCD as compared to the data [l5l |. 



Because of trigger bias, most of the contribution to di- 
hadron spectra comes from dijets produced close and tan- 
gential to the surface of the overlapped region, as shown 
in Fig. [3 However, there are still about 25% of the 
contribution coming from dijets near the center of the 
overlapped region. These jets truly punch through the 
medium and emerge after finite amount of energy loss. 
This is why dihadron spectra is slightly more sensitive to 
eo than the single spectra as was noted in Refs. fl^frTj ]. 
As one further increases the initial gluon density, the 
fraction of these punch-through jets will also diminish 
and the final dihadron spectra will be dominated by the 
tangential jets in the outer corona. Dihadron spectra at 
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FIG. 5: (color online). The same as Fig. [2] except for di- 
hadrons along the direction <j> = and n. 



RHIC will also lose its sensitivity to the initial gluon den- 
sity of the medium as shown in Fig.[3]Even if one includes 
the effect of transverse expansion of the bulk matter, the 
above results remain qualitatively the same [l7j . 

Also shown in Fig. [3] are the single and dihadron sup- 
pression factors at the LHC energy as a function of eo for 
fixed values of p^ ig = 20 GcV and pf? so = 10 GcV. Be- 
cause of the flattening of the overall spectra shape at the 
LHC energy, back-to-back dihadron spectra at a given 
Pt have more contribution from dijets with higher ini- 
tial energy. They are therefore less suppressed than at 
the RHIC energy and are more sensitive to the initial 
gluon density. The single hadron spectra on the other 
hand are increasingly dominated by surface emission and 
the suppression factor Raa becomes much smaller than 
Iaa for dihadrons. From a model estimate of the bulk 



hadron production at LHC [12j, the energy loss parame- 
ter is about eo ~ 5 GeV/fm in central Pb + Pb collisions. 
The dihadron suppression factor Iaa at around this value 
of eo is significantly larger than the single hadron sup- 
pression factor Raa and is more sensitive to the initial 
gluon density. 

In summary, we have used the NLO pQCD parton 
model with effective modified FF's due to radiative par- 
ton energy loss to study both single and dihadron spectra 
in high-energy heavy-ion collisions. We found that the 
surface emission dominates the single hadron production 
process within the range of initial gluon density in the 
most central Au + Au collisions at RHIC. However, there 
is still a significant fraction of dijets that are produced in 
the center of the dense medium and contribute to the fi- 
nal dihadron spectra after losing finite amount of energy. 
Therefore, dihadron spectra are more robust probes of 
the initial gluon density in the most central Au + Au col- 
lisions at RHIC while single hadron spectra become less 
sensitive. A simultaneous X 2 ~ n t to both the single and di- 
hadron spectra can be achieved within a narrow range of 
the energy loss parameters eo = 1.6 — 2.1 GeV/fm. If the 
initial gluon density at RHIC were to increase further, 
dihadron spectra at a fixed px would also be dominated 
by the surface emission and become insensitive to the 
initial gluon density. At LHC, however, the flattening of 
the initial jet production spectra leads to an increase in 
the dihadron suppression factor. Dihadrons will there- 
fore remain more robust probes within a large range of 
initial gluon density. 
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